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ABSTRACT 


A transient,  one-dimensional,  finite-difference  model  is  presented  that  predicts  the  coupled 
transfer  of  heat  and  moisture  in  a multilayer  wall  under  nonisothermal  conditions.  The  model 
can  predict  moisture  transfer  in  the  diffusion  through  the  capillary  flow  regimes.  It  has  a 
provision  to  account  for  convective  moisture  transfer  by  including  embedded  cavities  which  may 
be  coupled  to  indoor  and  outdoor  air. 

The  model  is  subsequently  used  to  predict  the  time-varying  average  moisture  content  in  the 
sheathing  and  siding  of  a wood  frame  wall  as  a function  of  time  of  year.  Results  are  shown  for 
a mild  winter  climate  (Atlanta,  GA),  an  intermediate  winter  climate  (Boston,  MA),  and  a cold 
winter  climate  (Madison,  WI).  The  indoor  temperature  is  maintained  at  21  °C,  and  separate 
computer  runs  are  carried  out  for  indoor  relative  humidities  of  35%  and  50%. 

The  effect  of  several  construction  parameters  on  the  winter  moisture  accumulation  are 
investigated.  The  parameters  include  the  interior  vapor  retarder  permeance,  sheathing 
permeance,  exterior  paint  permeance,  indoor  air  leakage,  and  the  amount  of  insulation. 


KEYWORDS:  capillary  transfer,  moisture  control  guideline,  moisture  management  in  walls, 
moisture  performance  of  walls,  moisture  transfer,  water  vapor  diffusion 


INTRODUCTION 

During  the  winter  season,  the  absolute  humidity,  or  moisture  content,  of  the  air  within  a 
residence  is  considerably  higher  than  that  of  the  outdoor  air.  As  a result,  moisture  permeates 
into  walls  by  way  of  diffusion  and  air  leakage  through  cracks  in  the  interior  surface.  This 
moisture  is  partially  adsorbed  and  accumulates  within  the  outer  material  layers  of  the  wall.  In 
response  to  the  seasonal  variation  in  outdoor  temperature,  Dutt  (1968)  observed  that  the  moisture 
content  in  outer  wall  layers  builds  up  during  cold  winter  periods  and  subsequently  decreases 
during  warm  summer  periods. 

Seasonal  moisture  cycling  of  the  outer  wall  layers  causes  wood-based  materials  to  alternately 
expand  and  contract.  Repeated  moisture  cycling  often  causes  warped  and  bowed  boards, 
delaminated  plywood,  pushed  out  nails,  and  the  separation  of  wood  members  from  the  structure. 
Moreover,  it  may  also  weaken  the  bond  between  exterior  paint  and  the  substrate  material, 
thereby  giving  rise  to  paint  failure.  In  isolated  cases,  high  moisture  contents  in  building 
materials  may  cause  fungus  degradation.  Moisture  problems  in  walls  have  been  documented  in 
field  surveys  carried  out  by  Tsongas  (1990),  Rose  (1986),  and  Merrill  and  TenWolde  (1989). 

Relative  to  the  above  discussion,  the  maximum  amount  of  sorbed  (or  bound)  moisture  that  can 
be  stored  in  a material  when  it  is  placed  in  an  environment  with  a relative  humidity  approaching 
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100%  is  denoted  "maximum  sorption".  Liquid  water  begins  to  appear  in  the  pore  structure  of 
the  material  when  its  moisture  content  approaches  and  exceeds  that  of  maximum  sorption.  The 
condition  existing  when  all  pore  structures  are  completely  filled  with  liquid  water  is  called 
"saturation".  For  conifer  softwoods  typically  used  in  building  construction,  the  moisture  content 
at  maximum  sorption  is  about  27%  of  the  dry  mass  of  wood  and  about  230%  at  saturation. 
Maximum  sorption  is  generally  regarded  as  the  maximum  amount  of  moisture  that  can  be  taken 
on  by  a material  without  degradation. 

Experts  often  disagree  on  the  cause  and  remedial  action  for  moisture  problems  in  residences. 
A contributing  factor  to  this  situation  is  that  general  analytical  models  have  not  been  available 
to  analyze  moisture  performance  of  building  components.  Mathematical  models  (e.g.,  Kohonen 
(1984),  Andersson  (1985),  Kiessl  (1983),  Oosterhout  and  Spolek  (1988),  and  Pedersen  (1990)) 
for  predicting  time-dependent  moisture  transfer  within  building  components  are  just  currently 
evolving.  Consequently,  experimentation  and  related  previous  experience  are  often  the  only 
proven  approaches  available  for  addressing  moisture  problems  while  minimizing  heat  loss. 
Experimentation  on  individual  components  is,  of  course,  costly  and  time-consuming.  Moreover, 
specific  results  cannot  be  readily  extended  to  different  constructions  and  indoor/outdoor  climatic 
conditions. 

Burch  and  Thomas,  et  al.  (1989)  previously  presented  a distributed-moisture-capacity,  finite- 
difference  model  that  predicted  the  time-dependent  moisture  diffusion  in  a multilayer  wall.  This 
model  was  one-dimensional  and  used  water  vapor  pressure  as  the  potential  for  moisture  transfer. 
The  model  was  experimentally  verified  for  a 3-layer  wall  tested  under  non-isothermal  conditions. 
This  previous  model  had  the  limitation  that  it  did  not  include  capillary  or  convective  transfer. 
In  the  present  paper,  the  earlier  finite-difference  model  is  extended  to  a more  general 
formulation  that  includes  these  effects. 


THEORY 

A composite  wall  comprised  of  N layers  in  series  that  stores  both  heat  and  moisture  is  analyzed 
(see  fig.  1).  Initially,  each  layer  has  an  arbitrary  moisture  content.  The  exterior  surfaces  of  the 
wall  are  then  exposed  to  ambient  environments  with  time-varying  temperatures  and  water  vapor 
pressures.  The  variation  in  moisture  content  and  temperature  within  each  of  the  layers  is  sought 
as  a function  of  time.  The  following  assumptions  are  used  in  the  analysis: 

• The  driving  forces  for  moisture  transfer  are  the  gradients  in  the  moisture  content  and 
temperature. 

• Heat  and  moisture  transfer  is  one-dimensional. 

• The  heat  transfer  properties  are  constant  (i.e.,  not  a function  of  temperature  or  moisture 
content). 
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The  sorption  isotherm  is  based  on  the  average  of  adsorption  and  desorption  data. 
Hysteresis  and  temperature  effects  on  the  sorption  isotherm  are  neglected. 

The  effects  of  temperature  and  hysteresis  on  the  relationship  between  suction  pressure  and 
moisture  content  are  neglected. 


Vapor  adsorption  at  a surface  releases  the  latent  heat  of  vaporization  and  vice  versa. 


Additional  assumptions  are  introduced  in  the  development  which  follows. 
Governing  Equations 


Within  each  layer  n of  the  wall  shown  in 
following  conservation  of  mass  equation^: 


dy 


figure  1,  moisture  transfer  is  governed  by  the 


dy 


Djiy.T) 


(1) 


The  selection  of  moisture  content  (7)  and  temperature  (T)  as  potentials  has  the  advantage  that 
the  same  mathematical  formulation  includes  both  diffusion  transfer  and  capillary  transfer,  as  will 
be  shown  later.  This  formulation  is  equivalent  (as  seen  later)  to  using  the  gradient  in  vapor 
pressure  as  the  moisture  transfer  potential  in  the  diffusion  regime  and  suction  pressure  in  the 
capillary  flow  regime,  with  a single  required  diffusivity. 


Heat  transfer  is  governed  by  the  conservation  of  energy  equation: 


by 


Kl.T) 


(2) 


Enthalpy  transport  by  moisture  movement  in  a wall  system  is  generally  small  and  therefore  is 
neglected  in  the  present  analysis.  However,  latent  transport  of  heat  is  included  in  the  boundary 
conditions  as  will  be  discussed  later.  The  term  (Q  +7C^)  includes  the  effect  of  energy  storage 
in  both  the  dry  material  and  accumulated  moisture.  The  effect  of  accumulated  moisture  can  be 
important  (e.g.,  the  specific  heat  of  wood  is  increased  by  69%  after  it  adsorbs  moisture  from 
air  at  a relative  humidity  of  100%). 


In  the  above  two  governing  equations,  strong  couplings  exist  between  heat  and  moisture  transfer. 
Both  the  diffusivity  for  the  moisture  gradient  (D^)  and  the  diffusivity  for  the  temperature 
gradient  (Dj)  are  strong  functions  of  moisture  content  and  temperature.  The  thermal 
conductivity  (k)  can  be  a function  of  moisture  content  and  temperature,  but  for  the  present 
analysis  it  is  assumed  to  be  constant. 


1 Symbols  are  defined  in  the  Nomenclature. 
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Indoor  Boundary  Conditions 


At  the  left  boundary  of  the  multilayer  wall  of  figure  1,  the  convective  heat  transfer  from  the 
indoor  air  plus  the  latent  heat  from  adsorbed  or  desorbed  water  vapor  is  equated  to  heat 
conduction  into  the  surface,  giving: 


- r)  + Tij  X = at  y ^ 0 (3) 

At  the  same  boundary,  the  moisture  transferred  through  an  air  film  and  paint  layer  is  equated 
to  moisture  transferred  into  the  surface,  or: 


M^Pi  - P)  = aty^O  (4) 

Here,  an  effective  conductance  (M^j)  defined  by: 

111 

= — + (5) 

Mei  Mj, 

has  been  introduced.  The  effect  of  a thin  vapor  retarder,  such  as  paint  or  a vinyl  covering,  is 
taken  into  account  as  a surface  conductance  (Mp^)  in  series  with  the  convective  mass  transfer 
coefficient  (M^j)  associated  with  the  air  film. 

In  equation  (4),  the  sorption  isotherm  function  (f),  defined  below,  is  used  as  a constitutive 
relation  to  evaluate  the  boundary  condition. 

7 =A4>,T)  ^M)  (6) 

A sorption  isotherm  is  illustrated  in  figure  2.  The  moisture  content  of  a porous  material  at  a 
relative  humidity  of  97%  may  be  measured  and  is  therefore  known.  Moreover,  when  a porous 
material  is  fully  saturated  with  liquid  water,  the  relative  humidity  within  the  pore  structure  is 
100%.  When  less  than  fully  saturated  with  water,  the  relative  humidity  is  less  than  100%  due 
to  the  presence  of  curved  menisci  within  the  pore  structure.  For  the  present  analysis,  a straight 
line  is  used  to  connect  the  moisture  contents  just  above  maximum  sorption  (<^  = 97%)  and  the 
state  of  saturation  (<^  = 100%). 

Outdoor  Boundary  Conditions 

Similarly,  at  the  right  boundary  shown  in  figure  1,  the  boundary  conditions  for  heat  and 
moisture  transfer  are: 
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h{T  -t',)*  til\  = 
and 

M.o[P  - Po)  = 

These  boundary  conditions  are  evaluated  in  a similar  fashion  as  those  at  the  indoor  surface  given 
in  equations  (3)  and  (4). 

Interface  between  Two  Storage  Lavers 

When  the  heat  transfer  at  the  interface  between  two  storage  layers  is  evaluated,  the  temperature 
is  assumed  to  be  continuous.  When  the  moisture  transfer  is  evaluated,  the  relative  humidity  is 
assumed  to  be  continuous,  and  the  sorption  isotherm  function  (f),  illustrated  in  figure  2,  provides 
the  constitutive  relationship,  or: 


■kH. 

dy 


at  y = L 


(7) 


- PdP>\ 


dy 


at  y = L 


(8) 


(9) 


where  the  subscripts  n and  n+1  refer  to  the  adjoining  regions.  The  relative  humidity  is 
continuous  at  an  interface,  but  the  moisture  content  will  generally  be  discontinuous. 

Non-Storage  Laver 


The  mathematical  model  has  the  provision  for  including  non-storage  layers  (e.g.,  an  air  space, 
a glass-fiber  insulation  cavity,  or  a vapor  retarder)  that  are  sandwiched  between  two  storage 
layers.  In  a non-storage  layer,  the  storage  of  heat  and  moisture  is  neglected,  and  the  transfers 
are  steady.  A non-storage  layer  may  be  convectively  coupled  to  indoor  and  outdoor  air. 


Consider  the  non-storage  layer  shown  in  figure  3.  At  the  interface  between  the  non-storage  layer 
and  the  adjacent  storage  layer  to  the  right,  the  heat  transferred  through  the  non-storage  layer  plus 
the  heat  gains  by  convective  exchange  with  the  indoor  and  outdoor  air  and  the  latent  heat  from 
adsorbed  or  desorbed  water  vapor  are  equated  to  the  heat  conduction  into  the  storage  layer,  or: 


s,n-l 


- T. 


s,n 


R 


vfpJT,  - r,,„)  - v.CpjT^  - 


dy 


at  y = Vs 


(10) 


For  the  moisture,  the  diffusion  transfer  through  the  non-storage  layer  plus  the  moisture  gains  due 
to  convective  exchange  with  the  indoor  and  outdoor  air  are  equated  to  the  moisture  transfer  into 
the  storage  layer,  or: 
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(11) 

r-i  dy  dT  . 

- ^y-ys 


The  sorption  isotherm  function  (f)  given  in  equation  (6)  is  used  as  a constitutive  relation  in 
evaluating  the  above  equation. 

Solution  Procedure 


0.622y,p,  0.622V^^ 

- fs,.)  - - Ps,)  - — 


Equations  (1-11)  were  recast  into  finite-difference  equations  using  a uniform  nodal  spacing 
within  each  layer.  An  implicit  solution  technique  with  coupling  between  the  two  conservation 
equations  was  used  to  solve  the  equations.  A FORTRAN  77  computer  program,  called  MOIST, 
with  a tridiagonal-matrix  solution  algorithm  was  prepared.  At  each  time  step,  the  calculation 
proceeds  by  first  solving  for  the  temperature  distribution.  Next,  a set  of  surface  moisture 
contents  are  calculated,  followed  by  a set  of  the  interior  moisture  contents.  Next,  a revised  set 
of  surface  moisture  contents  and  interior  moisture  contents  are  calculated  based  on  the  new 
values  of  the  diffusion  coefficients,  and  so  forth.  This  process  is  repeated  until  convergence  of 
the  moisture  contents  is  attained.  By  choosing  a sufficiently  small  time  step,  the  need  to  iterate 
between  the  temperature  and  moisture  solutions  was  eliminated. 

The  accuracy  of  the  numerical  solution  depends  on  both  the  nodal  spacing  and  the  time  step. 
Progressively  smaller  nodal  spacing  and  time  steps  were  used  to  ensure  that  calculated  results 
did  not  depend  significantly  on  the  space  or  time  increment. 

Later  in  the  paper,  program  MOIST  is  used  to  analyze  the  moisture  transfer  in  a wood  frame 
wall.  For  this  analysis,  a time  step  of  one  hour  is  used.  The  number  of  nodes  was  3 in  the 
gypsum  board,  5 in  the  sheathing,  and  5 in  the  wood  siding.  The  insulation  was  treated  as  an 
non-storage  layer.  When  the  computer  program  was  run  on  a Model  386  personal  computer, 
equipped  with  a math  co-processor  and  having  a 33  Mhz  clock  speed,  30  minutes  of  calculation 
time  was  required  to  simulate  one  year  of  real  time. 

MATERIAL  PROPERTIES 
Heat  Transfer  Properties 

The  thermal  conductivity,  density,  and  specific  heat  for  the  various  building  materials  used  in 
the  analysis  were  taken  from  ASHRAE  (1989)  and  are  summarized  in  table  1.  The  glass-fiber 
insulation  was  treated  as  a non-storage  layer  and  was  assumed  to  have  a thermal  resistance  of 
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1.9  m^*°C/W.  The  thermal  resistances  of  kraft  paper  and  the  paint  layers  were  neglected 
because  these  layers  are  very  thin. 

Diffusion  Regime 

In  the  regime  below  maximum  sorption,  free  liquid  water  is  not  present  in  the  pore  structure. 
Here  the  term  free  liquid  water  denotes  water  that  may  be  removed  by  body  forces  when  a 
material  is  spun  in  a centrifuge.  It  does  not  include  capillary  condensation  which  is  bound  in 
the  micropores  of  the  material.  Procedures  used  to  obtain  sorption  isotherms  and  vapor 
diffusivities  for  the  materials  are  discussed  below. 

Sorption  Isotherms.  The  sorption  isotherm  data  used  in  the  analysis  were  taken  from  Richards, 
et  al.(1992)  where  sorption  isotherm  data  were  fit  to  an  equation  of  the  form: 


(1  + a2(t>){l  - aj<l)) 

where  a^,  a2,  and  a3  are  empirical  constants  determined  by  a regression  analysis  of  the  measured 
data.  The  empirical  coefficients  are  summarized  in  table  2,  and  sorption  isotherm  curves  for 
the  materials  are  given  in  figure  4a. 

Moisture  Diffusivities.  The  water-vapor  diffusivities  for  the  materials  are  based  on  permeability 
measurements  carried  out  by  Burch,  et  al.  (1992).  For  each  material,  a series  of  permeability 
cup  measurements  were  carried  out.  In  the  measurements,  small  relative  humidity  differences 
were  imposed  across  the  specimens,  and  a functional  relationship  was  established  between 
permeability  and  the  average  relative  humidity  across  the  specimen.  Separate  measurements 
carried  out  at  24°  and  7°C  revealed  that  temperature  had  an  insignificant  effect  on  permeability. 
In  this  study,  the  permeability  data  were  fit  to  an  equation  of  the  form: 


/X  = exp(fli  + ^2^  + 

where  a^,  a2,  and  a3  are  empirical  constants  determined  from  a fit  of  the  measured  data.  The 
empirical  coefficients  are  summarized  in  table  3.  A plot  of  the  permeabilities  for  the  materials 
used  in  the  analysis  are  given  in  figure  4b. 

The  diffusivity  for  the  moisture  gradient  (D.^)  and  the  diffusivity  for  the  temperature  gradient 
(D-p)  were  calculated  by  the  relations: 


7 


and  D 


(14) 


mP^giT) 


Pd 


m) 

d(f> 


T ~ 


dP,g(T) 

dT 


Pd 


The  above  equations  may  be  derived  by  introducing  the  sorption  isotherm  function  and  applying 
the  chain  rule  to  Pick’s  steady-state  diffusion  equation  with  the  gradient  of  the  water- vapor 
pressure  as  the  driving-force  potential. 


The  glass-fiber  insulation,  the  19-mm-wide  air  space,  the  paint  layers,  and  the  kraft  paper  vapor 
retarder  were  treated  as  non-storage  layers.  Permeances  for  these  materials  are  based  on 
ASHRAE  (1989)  and  are  given  in  table  4. 


Capillary  Regime 


Capillary  transfer  occurs  when  a contiguous  path  of  liquid  exists  within  the  porous  material.  As 
the  moisture  content  of  a material  is  increased  above  maximum  sorption,  the  moisture  content 
at  which  a contiguous  path  of  liquid  first  exists  is  termed  "irreducible  saturation." 

Liquid  Diffusivitv.  In  the  capillary  flow  regime,  the  liquid  diffusivity  in  porous  materials  may 
be  predicted  by: 

SPc 

^ Pd^ 

This  equation  follows  by  applying  the  chain  rule  to  Darcy’s  Law  for  liquid  flow  through  a 
porous  media.  In  the  analysis,  the  density  of  water  (p^)  was  taken  to  be  1000  Kg/m^,  and  the 
viscosity  of  water  (v)  was  taken  to  be  7.25  X lO'"^  Pa^s.  Procedures  to  obtain  the  capillary 
pressure  (P^)  and  the  unsaturated  liquid  permeability  (k)  are  discussed  below. 


Capillary  Pressure.  The  Leverett  "j -function"  (Leverett  1941)  is  accepted  by  various  authors 
in  different  fields  as  a generalized  dimensionless  functional  form  that  may  be  used  to  correlate 
the  capillary  pressure  with  moisture  content  for  many  different  materials.  The  Leverett  j- 
function  is  defined  by: 


(16) 


where  is  the  liquid  permeability  of  the  porous  material  at  a saturated  moisture  content  and  a 
is  the  surface  tension  of  water,  taken  as  69.2  X 10'^  N/m.  Dry  porosity  (e^)  and  saturated  liquid 
permeability  values  (k^)  used  in  the  analysis  are  given  in  table  5. 


8 


The  Leverett  j -function  or  the  dimensionless  capillary  pressure  is  plotted  as  a function  of  the 
saturation  of  the  wetting  fluid  in  figure  5.  Here  the  saturation  of  the  wetting  fluid  (S)  is  defined 
as: 


S = 1 — (17) 

ys  - 

The  Leverett  j -function  was  curve  fitted  to  data  as  shown  in  figure  5. 

Unsaturated  Liquid  Permeability.  The  unsaturated  liquid  permeability  (k)  was  estimated  by 
the  linear  relation: 

K = kJS  (18) 


The  above  equation  is  based  on  the  modeling  work  of  Stanish,  et  al.  (1985).  Note  that  the 
un saturated  liquid  permeability  is  equal  to  zero  at  irreducible  saturation  and  equal  to  the 
saturated  liquid  permeability  (k^)  at  a fully  saturated  state. 

The  liquid  diffusivity  was  calculated  from  equation  (15)  using  the  above  procedures  to  estimate 
the  derivative  of  the  capillary  pressure  and  the  unsaturated  permeability.  The  authors  recognize 
that  this  procedure  provides  only  engineering  estimates  for  the  liquid  diffusivity  of  materials. 
However,  liquid  diffusivity  data  for  building  materials  is  seriously  lacking  in  the  literature. 
NIST  is  currently  measuring  liquid  diffusivities  for  building  materials.  As  these  measurement 
results  become  available,  they  will  be  incorporated  into  the  model  and  replace  the  approximate 
method  outlined  above. 

Within  the  capillary  flow  regime,  the  diffusivity  for  the  temperature  gradient  (D^)  was  calculated 
using  equation  (14). 

Transition  Regime 

When  the  moisture  content  of  a material  is  between  a state  of  maximum  sorption  and  irreducible 
saturation,  the  material  is  said  to  be  in  a "transition  regime."  At  and  below  the  irreducible 
saturation  moisture  content,  free  liquid  water  exists  but  not  in  a contiguous  path.  In  this  regime, 
the  capillary  attraction  between  discrete  liquid  particles  and  pores  is  so  strong  that  this  liquid 
cannot  be  separated  from  the  porous  material  by  ordinary  mechanical  means,  such  as 
centrifuging  or  applying  a pressure  gradient  across  the  material.  In  this  regime,  capillary 
transfer  vanishes  and  the  water  vapor  pressure  approaches  saturation. 

To  date,  different  researchers  do  not  agree  on  the  exact  moisture  transfer  mechanism,  especially 
in  an  isothermal  situation.  The  present  model  uses  moisture  concentration  and  temperature 
gradients  as  the  moisture  transfer  potential  in  this  regime.  The  diffusivity  for  the  moisture 
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gradient  (D^)  was  estimated  as  a function  of  moisture  content  by  joining  the  diffusivity  curves 
in  the  diffusion  and  capillary  regimes  with  a straight  line  (on  a logarithmic  scale)  as  shown  in 
figure  6a.  The  mathematical  relationship  for  the  diffusivity  in  the  transition  regime  is  then: 


D = D -Exp 

7 ms 


T Ims- 

Dir- 

_7,>*  - 7ms- 

Dms- 

(19) 


In  the  above  equation,  is  the  left  hand  limit  of  diffusivity  at  maximum  sorption  and 
is  the  right  hand  limit  of  liquid  diffusivity  at  irreducible  saturation.  Note  that  equations  (14), 
(15),  and  (19)  provide  a continuous  model  for  the  moisture  diffusivity  from  a dry  to  a saturated 
state. 

The  diffusivity  for  the  temperature  gradient  (Dj)  was  also  calculated  by  equation  (14)  and  is 
shown  in  figure  6b. 


DISCUSSION  OF  RESULTS 

The  above  mathematical  model  was  used  to  predict  the  seasonal  variation  in  the  moisture  content 
of  the  sheathing  and  wood  siding  of  the  wood  frame  wall  shown  in  figure  7.  This  wall  was 
comprised  of  13-mm-thick  gypsum  board  covered  with  latex  paint,  a 92-mm-thick  glass-fiber 
insulation,  13-mm-thick  fiber-board  sheathing,  and  13-mm-thick  wood  siding  covered  with  oil- 
base  paint. 

In  the  analysis,  the  indoor  temperature  was  maintained  at  21  °C.  The  outdoor  temperature, 
relative  humidity,  and  solar  radiation  were  taken  from  WYEC  hourly  weather  data  [Crow  1981] 
for  a mild  winter  climate  (Atlanta,  GA),  an  intermediate  winter  climate  (Boston,  MA),  and  a 
cold  winter  climate  (Madison,  WI).  The  heating  degree  days  for  these  cities  are  1,706°C ‘days 
for  Atlanta,  3207 °C  • days  for  Boston,  and  4228 °C  • days  for  Madison.  Weekly  average  outdoor 
temperatures  for  these  three  cities  are  given  in  figure  8.  Six  months  of  weather  data  were  used 
to  initialize  the  simulation  results  so  that  they  would  be  independent  of  assumed  initial  moisture 
content  and  temperature. 

Airtight  Wall  without  a Vapor  Retarder 

Winter  moisture  accumulation  in  an  airtight  (i.e.,  no  convection)  wood  frame  wall  without  a 
vapor  retarder  is  plotted  versus  time  of  year  in  figure  9a  for  an  indoor  relative  humidity  of  35  % 
and  in  figure  9b  for  an  indoor  relative  humidity  of  50%.  In  each  plot,  the  solid  line  depicts  the 
average  moisture  content  of  the  wood  siding,  while  the  broken  line  depicts  the  average  moisture 
content  in  the  fiber-board  sheathing.  In  the  figures,  maximum  sorption  for  the  wood  siding  is 
depicted  by  the  solid  horizontal  lines  and  that  for  the  sheathing  by  the  dashed  horizontal  line. 
As  indicated  earlier,  when  the  moisture  content  is  above  maximum  sorption,  free  liquid  water 
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exists  within  the  pores  of  a material,  and  serious  problems  are  expected  if  this  condition  exists 
over  an  extended  period  of  time. 

Comparing  the  pairs  of  curves  for  the  three  different  climates  in  figure  9a,  higher  moisture 
contents  are  seen  to  occur  in  colder  climates.  In  actual  houses,  the  effect  of  climate  will  be  less 
pronounced  because  the  indoor  relative  humidity  in  houses  tends  to  decrease  in  colder  climates 
as  a result  of  increased  moisture  losses  from  window  condensation  and  infiltration  of  drier  air. 
It  is  also  interesting  to  note  that  the  peak  moisture  content  occurs  about  two  months  after  the 
minimum  outdoor  temperature  shown  in  Fig  8. 

Comparing  Figs.  9a  and  9b,  indoor  relative  humidity  is  seen  to  be  a highly  significant  parameter 
affecting  winter  moisture  accumulation.  Increasing  the  indoor  relative  humidity  from  35%  to 
50%  substantially  increases  the  moisture  accumulation.  When  the  indoor  relative  humidity  was 
35%  (see  fig.  9a),  the  moisture  contents  of  the  sheathing  and  wood  siding  were  always  below 
maximum  sorption.  On  the  other  hand,  when  the  indoor  relative  humidity  was  increased  to  50% 
(see  fig.  9b),  the  moisture  contents  approached  maximum  sorption  in  Boston  and  rose 
considerably  above  maximum  sorption  for  4-6  months  in  Madison. 

These  large  seasonal  variations  in  moisture  content  will  cause  wood-based  products  to  undergo 
expansion  and  contraction  cycles  which  contribute  to  material  degradation. 

The  results  in  figure  9 and  with  the  same  type  of  results  following  show  that  peak  moisture 
accumulation  in  various  materials  generally  do  not  exceed  maximum  sorption  by  much  and,  in 
fact,  never  reach  the  capillary  flow  regime.  This  situation  means  that  the  overall  results  are 
relatively  insensitive  to  the  significant  uncertainties  in  material  properties  above  maximum 
sorption. 

Airtight  Wall  with  a Vapor  Retarder 

A vapor  retarder,  placed  between  the  gypsum  board  and  the  insulation,  was  added  to  the  wall 
used  in  the  two  previous  computer  simulations.  The  vapor  retarder  consisted  of  a kraft  paper 
having  a permeance  of  1 .7  X 10'^^  kg/s  • nr  • Pa.  The  simulation  results  with  this  vapor  retarder 
are  given  in  figure  10  for  the  same  two  indoor  relative  humidity  conditions. 

Comparing  Figs.  9 and  10,  the  vapor  retarder  is  seen  to  provide  a significant  reduction  in  the 
seasonal  cycles  in  moisture  content.  Note  that  the  moisture  content  in  both  the  sheathing  and 
wood  siding  is  always  maintained  below  12%  even  for  the  higher  indoor  relative  humidity.  This 
reduces  the  seasonal  expansion  and  contraction  cycles  in  wood-based  products,  thereby  reducing 
material  degredation  due  to  moisture  cycling.  These  results  indicate  that  a vapor  retarder  is  an 
effective  means  of  controlling  moisture  accumulation  in  a wood  frame  wall  exposed  to  cold 
climates. 
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Wall  with  a Vapor  Retarder  and  Indoor  Air  Leakage 


The  mathematical  model  was  next  used  to  investigate  the  effect  of  indoor  air  leakage  on  the 
moisture  accumulation  within  the  wood  frame  wall  with  a vapor  retarder  (see  fig.  11).  In  the 
analysis,  the  insulation  cavity  was  ventilated  with  indoor  air  at  an  assumed  rate  of  one  volume 
change/h.  Convective  air  exchange  occurs  due  to  air  leakage  paths  associated  with  electrical 
outlet  boxes,  baseboard  cracks,  and  other  inherent  penetrations. 

The  moisture  contents  of  the  wood  siding  and  fiber-board  sheathing  are  plotted  versus  time  of 
year  in  figure  12a  for  an  indoor  relative  humidity  of  35%  and  in  figure  12b  for  an  indoor 
relative  humidity  of  50%.  For  an  indoor  relative  humidity  of  50%,  the  moisture  contents  are 
seen  to  approach  maximum  sorption  in  Madison.  It  should  be  pointed  out  that  this  analysis  is 
one-dimensional  and  does  not  include  the  effect  of  higher  sheathing  moisture  contents  directly 
in  line  with  impinging  air  leakage.  The  effect  of  convective  transport  is  averaged  over  the  entire 
sheathing  surface. 

The  results  of  figure  12  suggest  that  indoor  air  leakage  is  also  an  important  moisture  transfer 
mechanism.  Comparing  the  results  of  Figs.  10  and  12,  indoor  air  leakage  is  seen  to  substantially 
increase  the  moisture  content  in  a wall  with  a typical  vapor  retarder.  These  results  indicate  that 
in  order  to  reduce  moisture  accumulation  in  a wood  frame  wall  exposed  to  cold  climate,  it  is 
also  necessary  to  seal  interior  air  leakage  paths. 

It  should  be  pointed  out  that  the  above  analysis  is  very  dependent  on  an  assumed  cavity 
ventilation  rate  of  one  volume  change/h.  This  value  was  not  based  on  measured  data,  but  rather 
an  engineering  estimate  that  seemed  reasonable  to  the  authors. 

Effect  of  Other  Construction  Parameters 


In  this  section,  a sensitivity  analysis,  using  the  wood  frame  wall  of  figure  7,  is  carried  out  to 
investigate  the  effect  of  orientation,  exterior  paint  permeance,  exterior  sheathing  permeance,  and 
the  amount  of  insulation.  In  the  analysis,  the  wall  is  exposed  to  an  indoor  relative  humidity  of 
35%  and  the  outdoor  weather  of  Madison,  WI. 

Orientation.  The  moisture  contents  of  the  sheathing  and  siding  of  a north-facing  and  south- 
facing wall  are  compared  in  figure  13a.  The  higher  solar  radiation  for  the  south-facing  wall 
causes  its  moisture  contents  to  be  lower  than  those  for  an  identical  north-facing  wall.  The  peak 
wood  moisture  content  in  the  south-facing  wall  is  0.04  lower  than  for  the  identical  north-facing 
wall. 

Exterior  Paint  Permeance.  In  figure  13b.  the  moisture  accumulation  within  the  wood  frame 
wall  with  exterior  oil-base  paint  is  compared  to  the  same  wall  with  a considerably  more 
permeable  latex  paint.  The  permeance  of  the  latex  paint  was  5.7  X 10‘^°  kg/s  • m^  • Pa,  and  that 
for  the  oil-base  paint  was  1.1  X 10'^^  kg/s ‘m"*  Pa. 
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The  results  indicate  that  the  use  of  exterior  latex  paint,  as  opposed  to  oil-base  paint,  enhances 
the  escape  of  moisture  during  warm  drying  periods,  and  thereby  maintains  somewhat  lower 
moisture  contents  in  the  wood  siding.  The  peak  wood  moisture  content  in  the  wall  painted  with 
latex  paint  is  seen  to  be  0.03  lower  than  the  same  wall  painted  with  oil-base  paint. 

Sheathing  Permeance.  The  moisture  accumulation  for  a wood  frame  wall  with  fiber-board 
sheathing  is  compared  to  that  for  the  same  wall  with  the  considerably  less  permeable  plywood 
sheathing.  The  permeance  of  the  fiber-board  sheathing  was  about  150  times  that  of  the  plywood. 
The  results  are  given  in  figure  13c. 

The  results  indicate  that  the  use  of  less  permeable  plywood  sheathing  reduced  the  ingress  of 
moisture  into  the  wood  siding.  The  peak  moisture  content  in  the  wood  siding  for  the  wall  with 
plywood  sheathing  was  0.05  lower  than  that  for  the  same  wall  with  fiber-board  sheathing.  Note 
that  the  moisture  content  of  the  two  sheathing  materials  do  not  differ  much. 

Amount  of  Insulation.  In  figure  13d,  the  moisture  accumulation  in  a wood  frame  wall  without 
insulation  is  compared  to  the  same  wall  with  insulation.  When  the  thermal  insulation  is  in  the 
wood  frame  wall,  most  of  the  inside- to-outside  temperature  difference  occurs  across  the 
insulation,  and  a small  temperature  difference  occurs  between  the  sheathing  and  wood  siding. 
This  causes  the  moisture  content  of  the  sheathing  and  siding  not  to  be  much  different.  On  the 
other  hand,  when  the  thermal  insulation  is  absent  from  the  wall,  a much  larger  temperature 
difference  occurs  between  the  sheathing  and  wood  siding.  This  causes  moisture  to  move  from 
the  warmer  sheathing  outwards  into  the  colder  wood  siding.  This  result  indicates  that  when 
uninsulated  wood  frame  walls  are  retrofitted  with  thermal  insulation  in  cold  climates,  the 
sheathing  materials  will  tend  to  have  increased  moisture  content. 

Using  Outdoor  Ventilation  as  a Moisture  Management  Technique 

Another  moisture  management  technique  for  wood  frame  walls  is  the  use  of  outdoor  ventilation 
to  remove  accumulated  moisture  instead  of  using  an  interior  vapor  retarder  and  sealing  interior 
air  leakage  paths.  Here  a 19-mm  naturally-ventilated  cavity  is  formed  between  the  sheathing  and 
wood  siding  as  shown  in  figure  14.  In  the  analysis,  the  cavity  ventilation  rate  was  assumed  to 
be  6 volume  changes/h. 

The  results  of  the  analysis  are  given  in  figure  15a  for  an  indoor  relative  humidity  of  35%  and 
in  figure  15b  for  an  indoor  relative  humidity  of  50%.  In  each  figure,  the  following  three  cases 
are  compared:  the  original  wood  frame  wall  shown  in  figure  7,  the  wood  frame  wall  with  an 
unventilated  air  cavity,  and  the  wood  frame  wall  with  the  ventilated  air  cavity  shown  in  figure 
14. 

These  results  show  that,  when  the  sheathing  and  siding  are  in  contact,  the  moisture  content  of 
the  sheathing  and  wood  siding  are  not  much  different.  When  an  unventilated  air  space  is  placed 
between  the  sheathing  and  wood  siding,  the  moisture  content  of  the  wood  siding  increases  and 
the  moisture  content  of  the  sheathing  decreases.  This  effect  is  caused  primarily  by  the  thermal 
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resistance  of  the  air  space  altering  the  temperature  distribution  within  the  construction  (i.e.,  the 
wood  siding  becomes  colder  and  the  sheathing  becomes  warmer).  The  larger  temperature 
difference  between  the  wood  siding  and  sheathing  causes  accumulated  moisture  to  redistribute 
towards  the  colder  wood  siding.  When  the  air  space  is  ventilated,  the  moisture  content  of  both 
the  sheathing  and  the  wood  siding  decrease. 

Comparing  Figs.  15  and  10,  it  is  seen  that  the  use  of  outdoor  ventilation  is  considerably  less 
effective  in  reducing  seasonal  fluctuations  in  moisture  content  compared  to  the  practice  of 
providing  a vapor  retarder  and  sealing  air  leakage  paths  at  the  interior  surface. 


SUMMARY  AND  CONCLUSIONS 

A mathematical  model  that  predicts  the  combined  transfer  of  heat  and  moisture  in  walls  was 
presented.  The  model  is  one-dimensional  and  includes  diffusion,  capillary,  and  convective 
transfer.  The  model  was  used  to  analyze  the  winter  moisture  accumulation  in  a wood  frame  wall 
exposed  to  a mild  winter  climate  (Atlanta,  GA),  an  intermediate  winter  climate  (Boston,  MA), 
and  a cold  winter  climate  (Madison,  WI). 

Seasonal  changes  in  outdoor  temperature  were  found  to  produce  seasonal  variations  in  the 
moisture  content  of  the  outer  layers  of  a wood  frame  wall.  The  most  important  parameters 
affecting  the  amount  of  moisture  accumulation  during  the  winter  were  the  indoor  relative 
humidity  and  the  outdoor  climate.  Indoor  relative  humidity  was  observed  to  be  more  important 
than  outdoor  climate.  The  amount  of  moisture  accumulation  was  larger  in  colder  climates. 
Other  parameters  found  to  have  a less  important  effect  on  the  moisture  accumulation  were  wall 
orientation,  exterior  paint  permeance,  and  sheathing  permeance.  The  addition  of  insulation  to 
an  uninsulated  wood  frame  wall  was  found  to  significantly  increase  the  moisture  content  of  fiber- 
board  sheathing. 

Both  diffusion  and  indoor  air  leakage  were  shown  to  be  important  mechanisms  for  transporting 
moisture.  In  walls  without  an  interior  vapor  retarder,  moisture  diffusion  was  shown  to 
accumulate  sufficient  moisture  to  produce  free  liquid  water  in  the  sheathing  and  wood  siding  of 
an  airtight  wall  exposed  to  an  indoor  relative  humidity  of  50%  and  outdoor  climate  of  Madison, 
WI.  The  leakage  of  indoor  air  into  the  insulation  cavity  of  a wood  frame  wall  with  an  interior 
vapor  retarder  was  observed  to  accumulate  moisture  that  approached  a state  of  free  liquid  water. 

The  inclusion  of  a vapor  retarder  and  sealing  air  leakage  paths  at  the  interior  wall  surface  were 
found  to  decrease  significantly  seasonal  variations  in  moisture  contents,  thereby  preventing  large 
expansion  and  contraction  cycles  in  wood-based  products  that  are  believed  to  degrade  the 
construction.  The  practice  of  using  outdoor  ventilation  in  a cavity  between  the  sheathing  and 
siding  was  found  to  be  less  effective  than  the  practice  of  including  a vapor  retarder  and  sealing 
interior  air  leakage  paths. 
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NOMENCLATURE 

Symbol  Units 

C J/kg-°C 

m^/s 

Dt  kg^^,-m-/kgj-°C-s 
f kgjkgd 

h W/m^-X 

j 

k W/m  • X 

L m 

M kg^s  • m“  • Pa 

h " kg!  VC?  • s 

P Pa 

R m-*°C/W 

S 

t s 


y m 

V VC? Is  per  m^ 

7 kgjkgd 

6 

2 

K m 

X J/kg^ 

/X  kg^s*m^*Pa 

V Pa*s 

p kg/m^ 

<j  N/m 

4> 


Definition 


Constants  in  sorption  isotherm  and  permeability  function 
(n=l,2,  and  3) 

Specific  heat 

Diffusivity  for  moisture  gradient 
Diffusivity  for  temperature  gradient 
Sorption  isotherm  function 
Surface  heat  transfer  coefficient 
Dimensionless  capillary  pressure 
Thermal  conductivity  of  porous  material 
Thickness  of  wall 

Permeance  or  moisture  conductance 

Moisture  mass  flux 

Pressure 

Thermal  resistance 

Saturation  of  wetting  fluid  (see  equation  17) 

Time 

Temperature 

Sol-air  temperature 

Distance  from  inside  surface  of  wall 

Volumetric  air  flow  rate  per  unit  area 
Moisture  content  on  dry  basis 
Porosity 

Un saturated  liquid  permeability 
Latent  heat  of  vaporization 
Water-vapor  permeability 
Viscosity  of  water 
Density 

Surface  tension  of  water 
Relative  humidity 
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Subscripts  Refer  to: 


a = Atmospheric  or  air  property 
c = Capillary 

d = Dry  property 

e = Effective  property 
f = Air  film 

g = Saturated  state 

i = Indoor  or  inside  surface  property 

ir  = Irreducible  saturation 
ms  = Maximum  sorption 
n = Storage  layer  index 
o = Outdoor  property 
p = Paint  property 

s = Surface  value  or  liquid  saturated  state 

T = Temperature  gradient 
V = Vapor  property 

w = Moist  or  water  property 
7 = Moisture  content  gradient 
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TABLE  1 

Heat  Transfer  Properties 

Material 

Thermal 

Density 

Specific  Heat 

Conductivity 

kg/m^ 

J/kg-”C 

Plywood 

0.115 

545 

1214 

Wood  Siding 

0.118 

365 

1633 

Gypsum  Board 

0.160 

670 

1098 

Fiber-Board  Sh. 

0.0663 

352 

1298 

TABLE  2 

Empirical  Constants  for  Sorption  Isotherm  Function 

Material 

aj 

^2 

^3 

Plywood 

0.3441 

6.177 

0.8283 

Wood  Siding 

0.1936 

2.095 

0.7687 

Gypsum  Board 

0.02465 

9.0750 

0.93540 

Fiber-Board  Sh. 

2.7136 

65.267 

0.8684 

TABLE  3 

Empirical  Constants  for  Permeability  Function 

Material 

as 

Plywood 

-30.10 

3.2963 

2.4391 

Wood  Siding 

-28.68 

-0.9198 

4.5776 

Gypsum  Board 

-23.470 

-1.4799 

1.0816 

Fiber-Board  Sh. 

-24.09 

0.1836 

-0.3919 
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TABLE  4 

Permeances  for  Non-Storage  Layers 

Material 

10'^®  kg^/s*m^‘  Pa 

Latex  Paint 

5.7 

Oil-Base  Paint 

1.1 

92  mm  Glass  Fiber  Ins. 

18. 

Kraft  Paper  (Asphalt  Impreg.) 

0.17 

19  mm  Air  Space 

69. 

TABLE  5 

Dry  Porosity  and  Saturated  Liquid  Permeability  Values  Used  in  Analysis 

Material 

. -'O' 

m^ 

Plywood 

0.636 

1.4  X 10'“ 

Wood  Siding 

0.756 

2.8  X 10-*’ 

Fiber-Board  Sh. 

0.765 

3.9  X 10-*5 
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MOISTURE  CONTENT  (Y) 


•0’  '”60 


Fig.  1 . Wall  consisting  of  N layers 


Fig.  2.  The  sorption  isotherm 
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Fig  3.  A non-storage  layer 
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RELATIVE  HUMIDITY 

b.  Permeabilities 

Fig.  4.  Diffusion  properties  of  materials  used  in  the  analysis 
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DIMENSIONLESS  CAPILLARY  PRESSURE  (J) 
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0.2  0.4  0.6  0.8  1.0 


SATURATION  OF  WETTING  FLUID  (S) 


Fig.  5.  Dimensionless  capillary  pressure  versus  saturation 

of  wetting  fluid  for  unconsolidated  sands  [Collins  (1961)] 
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, m^/s 


a.  Diffusivity  for  the  moisture  content  gradient 


Fig.  6.  Moisture  diffusivity  for  white  pine 
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13  mm  fiber  board  sheathing 

13  mm  wood  siding 


Oil-base  paint 


Fig.  7.  Wood  frame  wall  used  in  the  analysis 


TIME  OF  YEAR 


Fig.  8.  Weekly  average  outdoor  temperature  for  the 
three  outdoor  climates 
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TIME  OF  YEAR 

a.  indoor  relative  humidity,  (|)  = 35% 


TIME  OF  YEAR 

b.  Indoor  relative  humidity,  (|)  = 50% 

Fig.  9.  Moisture  content  plotted  versus  time  of  year  for  the  wood 
frame  wall  without  a vapor  retarder  and  without  indoor 
air  leakage 
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MOISTURE  CONTENT  MOISTURE  CONTENT 


TIME  OF  YEAR 

a.  Indoor  relative  humidity,  ([)  = 35% 


TIME  OF  YEAR 

b.  Indoor  relative  humidity,  <|)  = 50% 


Fig.  1 0.  Moisture  content  plotted  versus  time  of  year  for  the  wood 
frame  wall  with  a vapor  retarder  and  without  indoor 
air  leakage 
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MOISTURE  CONTENT 


Fig.  1 1 . Illustration  of  indoor  air  leakage  into  insulation  cavity 
of  wood  frame  wall 


TIME  OF  YEAR  TIME  OF  YEAR 


a.  Indoor  relative  humidity,  (j)  = 35%  b.  Indoor  relative  humidity,  (|)  = 50% 


Fig.  1 2.  Moisture  content  plotted  versus  time  of  year  for  the  wood  frame 
wall  with  a vapor  retarder  and  with  indoor  air  leakage 
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MOISTURE  CONTENT  MOISTURE  CONTENT 
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a.  Effect  of  orientation  b.  Effect  of  exterior  paint  permeance 


TIME  OF  YEAR 


TIME  OF  YEAR 


c.  Effect  of  sheathing  permeance  d.  Effect  of  cavity  insulation 


Fig.  1 3.  Moisture  content  plotted  versus  time  of  year  for  the  wood  frame 
wall  without  a vapor  retarder  and  without  indoor  air  leakage 
(Madison,  Wl.) 
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MOISTURE  CONTENT 


Fig.  14.  Illustration  of  wood  frame  cavity  wall  with  an  exterior 
cavity  ventilated  with  outdoor  air 


a.  Indoor  relative  humidity,  (j)  = 35%  b.  Indoor  relative  humidity,  (j)  = 50% 


Fig.  15.  Comparison  of  moisture  content  variation  of  wood  frame  wall 
without  outdoor  ventilation  (Fig.  7)  versus  that  of  the  same 
wood  frame  wall  with  outdoor  ventilation  (Fig.  14)  in  Madison  Wl. 
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